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Abstract We have obtained the spectrum of a middle- 
aged PSR B0656+14 in the 4300-9000A range with the 
ES0/VLT/F0RS2. Preliminary results show that at 
4600-7000A the spectrum is almost featureless and flat 
with a spectral index 0^:^-0.2 that undergoes a change to 
a positive value at longer wavelengths. Combining with 
available multiwavelength data suggests two wide, red 
and blue, flux depressions whose frequency ratio is about 
2 and which could be the 1st and 2nd harmonics of elec- 
tron/positron cyclotron absorption formed at magnetic 
fields ~10 8 G in upper magnetosphere of the pulsar. 
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1 Introduction 

Multiwavelength observations of radio pulsars are an 
important tool for the study of not yet clearly un- 
derstood radiative mechanisms and spectral evolution 
' of rotationpowered isolated neutron stars (NSs). Opti- 
cal observations are an essential part of these studies. 
Among a dozen optically identified NSs, only seven pul- 
sars have parallax based distances and, therefore, min- 
imal uncertainties (<10%) in luminosities. Their op- 
tical spectra are shown in Fig. 1. A reliable optical 
spectrum has been obtained only for th e young and 
bright Crab pulsar l)Sollerman et al.lf2000(l . Other pul- 
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Fig. 1 Optical spectra of seven pulsars of different character- 
istic age indicated in the left-bottom corner of each panel in 
years. The youngest Crab pulsar is at the top and the oldest 
PSR B0950-08 is at the bottom. For PSR B0656+14 dashed 
lines show the low energy extensions of the blackbody (BB, 
T=0.84MK) and power la w (PL, a = 0.45) X-ray spectral 
components and their sum (JKoptsevich ct al. 20(21). 



sars are fainter and are represented mainly by broad- 
band photome tric points. Published optical spectra o f 
PSR B0540-69 ijHill et al.ll997tlSeranmovich et al.l200^ 
are strongly contaminat ed by a bright pulsar nebula 
l)Seranmoyich et aT1l2004lh wh ile a tentative spectrum of 
Geminga~ i|Martni^t^dTir998j) is much n o isier t han avail- 
able photometric nuxesT lMignam et alJ l)2004|) reported 
on the Vela pulsar spectral observations but these data 
are not published yet. 
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. The correlation coefficient is 0.97. 
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Fig. 2 From top to bottom: Evolution of the spindown, ra- 
dio, optical, X-ray, and 7-ray luminosities and respective ef- 
ficiencies demonstrated by the 7 optical pulsars from Fig. 1. 



In Figure 2 we show the evolution of luminosity L 
and radiation efficiency r\ = L/L s( j (L s d is spindown lu- 
minosity) demonstr ated by these seven pu lsars using the 
data of Table 1 from lZharikov et all l)2006[) . We note sig- 
nificantly non-monotonic dependencies of r/opt and r/x 
versus pulsar age with a pronounced minimum at the 
beginning of the middle-age epoch (~ 10 yr) and com- 
parably higher efficiencies of younger and older pulsars. 

Owing to its relative brightness and proximity, the 
middle-aged PSR B0656+14 is one of isolated NSs most 
intensively studied in different wave l ength s. It was dis- 
covered in radio bv lManchester et alJ l)l978|) . then identi- 
fied in X-rays with Einstein, and observed in details with 
ROSAT, ASCA , Chandra and XMM (for references see 
IShibanov et all 1120051 1200(^ 1. The X-ray emission can 
be described as a combination of thermal radiation from 
the entire surface of a cooling NS and from hotter polar 
caps heated by relativistic particles of magnetospheric 
origin. An excess over the hot thermal component at 
energies >2 keV was interpreted a s nonthermal radia- 
tion f rom the pulsar magnetosphere l)Greiveldinger et all 
1996). The pulsar ha s been also marginally detect ed in 
7-rays (>50 MeV) bv iRamanamurthv et all l(l99rJ ). 

In the op ti cal P SR B0656+14 was identified by 
ICaraveo et all l|l994h with the ESO/NTT telescopes 
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Fig. 3 Relationship between the B-band optical and 2-10 
keV X-ray efficiencies for the same 7 pulsars as in Fig. 1,2. 



in the V band. It was then studied in UV with the 
HST/F OC in the F130LP F430W, F342W and F195W 
bands l|Pavlov et all 119971) and with t he H ST/WFPC 
in the F555W band l|Mignani et all Il997|) . Detailed 
photometr ic studies in the optical-NIR were per - 
formed bylKurt et all l ll998|). iKootsevich et aH (l200lfl . 
iKomarova et all l)2003|) and lshibanov et all 1)20061) . The 
studies showed that the bulk of the optical radiation is 
of nonthermal origin. This was confirmed by the detec- 
tion of coherent optical pulsations with the radio pulsar 
period in the B band SJiga^^talJ_M9J3) , in a wide 
400-600 nm passband ijKern et all l2003[) and in NUV 
( Shib anov et all l2005|h The pulse profiles are rather 
sharp with a high pulse fraction as expected for nonther- 
mal emission mechanisms. The phase integrated multi- 
wavel ength spectrum (Fig. 4) shows ijKoptsevich et all 
120011 IShibanov etHI [2006) that the NIR-optical-UV 
spectral energy distribution is, in a first approach, com- 
patible with the low energy extension of the sum of the 
X-ray thermal blackbody (BB) spectral component from 
the whole NS surface and the power law (PL) compo- 
nent dominating in the high energy tail. The BB ex- 
tension does not contribute at longer wavelengths where 
the optical-NIR fluxes are in a good agreement with the 
PL alone (Fig. 1). This indicates a common origin of the 
nonthermal optical and X-ray emission, which is strongly 
supported by a good coincidence in phase and shape 
of the pulse profiles in the optical and the X-ray tail 
l|Shibanov et al.ll2005j) . The same origin of the nonther- 
mal optical and X-ray photons is likely to be a general 
property for other pulsars detected in both ranges, as 
follows from a strong corr elation between respective ef - 
ficiencies (Fig. 3) found bv lZharikov et all l)2004l2006[) . 

In the optical emission of PSR B0656+14 there is an 
apparent, ~(3-5)<r, flux excess over the PL "continuum" 
at Log(i^) ~14.7 (Fig. 1). This could reveal an additional, 
3rd, spectral component to the BB+PL discussed above. 
Here we present first results of the optical spectroscopy 
of the pulsar partially motivated by more detailed studies 
of the excess. In a broader sense, these results also allow 
us to consider, for the first time, optical properties of a 
middle-aged ordinary pulsar at the spectroscopic level, 
as it has been achieved so far only for the Crab pulsar. 
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Fig. 4 Left: Unabsorbed multiwavelength spectrum of PSR B0656+14 from the radio through 7-rays. The box marks the 
range zoomed in the right panel. Right: The spectral and photometric fluxes of PSR B0656+14 in the NIR-optical-NUV range 
obtained with different telescopes and instruments, as notified in the plot. Black dashed lines show low energy extensions of 
the soft blackbody (BB) and power law (PL) spectral fits and their sum (BB+PL) obtained in X-rays. The black solid lines 
show the same but with the PL normalization shifted up by a factor of 1.4 (PL1) to fit the upper edge of its la error bar 
shown at the left side of the plot. These are possibly a better match for the optical and NUV spectra than the dashed lines. 
The symbol © marks the Earth atmospheric absorption band near 7600A. 



2 Observations and data analysis 

The spectrum of PSR B0656+14 was obtained on 
November-December 2004 and February 2005 during 
several observational runs of the ESO program 074. D- 
0512A using the VLT/UT1 telescope in a service mode. 
The FORS2 1 instrument was used in a long slit spectro- 
scopic setup with the grating GRIS_300V and the filter 
GG435, which cover the wavelength interval of about 
4300-9600A and provide a medium spectral resolution of 
3.35A/pixel. The slit width was 1" and its position angle 
was selected in a such way as to obtain also spectra of 
several nearby stars for a sure pulsar astrometric refer- 
encing and wavelength/flux calibration. Eighteen 1400 s 
science spectroscopic exposures were taken with a total 
exposure time of 25200 s at a mean seeing of 0.6". Stan- 
dard reference frames (biases, darks, flatfields, lamps) 
were obtained in each observational run, while the slit 
and slitless observations of spectrophotometric standards 
(FeigellO, LTT3218 and LTT1788) for the flux calibra- 
tion were carried out in separate runs on the same nights. 

A combination of the MIDAS and IRAF packages 
was used for standard CCD data reduction, cosmic-ray 
track removing, spectra extraction, and subsequent data 
analysis. A faint, R=24.65, pulsar is at a limit of spectro- 
scopic capability of the VLT. Nevertheless, excellent see- 
ing conditions allowed us to resolve its spectrum even at 
each individual exposure, albeit with a low signal to noise 
ratio S/N. These exposures were co-added. The spec- 
trum was then extracted with a 3 pixel wide extraction 
slit (0.2" /pix) centered on the pulsar. The backgrounds 
were extracted with a 6 pixel wide slit centered above and 
below the center of the pulsar spectrum. The correction 

1 www.eso.org/instruments/fors 



factor for the PSF and sensitivity function were obtained 
from the FeigellO standard observations. The S/N of the 
resulting spectrum was about 4 (per pixel) in the 4450- 
5500A range and declined to ~ 1 near/above 8000A, due 
to higher sky backgrounds and a drop in sensitivity to- 
wards longer wavelengths. We binned the spectral flux 
in 20 pixel bins (67A) to get S/N near/above 15 and 4, 
respectively, making the flux accuracy to be comparable 
with that of available photometric data. 



3 Results and discussion 

The binned and dereddened with Ay=0.093 spectrum of 
the pulsar is shown in Figure 4 (red curve in the right 
panel) . We show also a vailable multiwavelength data (see 
IShibanov et alJ l)2006|) for the data and Ay description). 
The spectrum is in a good agreement with the broadband 
VRI fluxes, while it is somewhat higher than the B band 
flux. It does not show any strong nebular emission lines 
that could be responsible for the apparent VR excess, 
mentioned above as a 3rd component, while the presence 
of weak features can not be completely ruled out. The 
continuum of the dereddened spectrum in 4600-7000A 
range has a power law shape oc j,-°- 2 ± - 2 (green line) 
in agreement with the nonthermal nature of the bulk 
of the optical emission. Its slope is close to the value 
expected in this range from the BB+PL extended from 
X-rays. The slope likely undergoes a change to a positive 
value between the R and I bands, as follows also from the 
photometric data. 

Within uncertainties the bluest end of the optic al and 
the reddest end of the NUV l|Shibanov et al.ll20 05) spec- 
tra are compatible with each other, suggesting a smooth 
connection of both. However, the F430W photometric 
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flux in the gap between them drops below this connec- 
tion with a significance of about 2a of the flux confidence 
level. Unless this is a result of some unknown systemat- 
ics in calibration, it suggests a spectral dip in the pulsar 
emission centered near -4300A (-14.83 in Log(z/)). 

The optical and NUV spectra and two NIR pho- 
tometric points, F160W and F187W, almost perfectly 
match the BB+PL extension from X-rays, if the PL nor- 
malization is taken to be a factor of 1.4 higher (solid 
lines) than its best X-ray fit value (dashed lines). The 
change is within la uncertainty of the fit. Consider- 
ing the solid line version as a new optical "continuum 
level" we find an additional and more significant flux 
depression in the red part of the spectrum overlapping 
the I and F110W bands and centered near 9000-10000A 
(Log(i/)~14.5). 

Additional spectral studies are necessary to confirm 
the suggested "red" and "blue" features and to mea- 
sure their shapes and wavelengths more accurately. Nev- 
ertheless, if they are real, an approximate bluc-to-rcd 
frequency ratio is ~2. This indicates that they can be 
the 1st and 2nd harmonics of an electron/positron cy- 
clotron absorption formed in the upper magnetosphcrc 
of the pulsar at an effective altitude where the magnetic 
field B~10 8 G. This is ~360 km, assuming a dipole NS 
field with a surface value of 4.66xl0 12 G, as derived 
from spindown measurements. The absorbing e have to 
be cooled enough and provide a sufficient optical depth 
above a source of the nonthermal continuum. The source 
altitude is, therefore, <360 km and much below the light 
cylinder radius of ~18xl0 3 km, likely suggesting its po- 
lar cap origin. The features are broad, as expected from 
the magnetospheric field inhomogencity. Tentative (<2a) 
absorption features in the NUV s pectrum at Log(^) ~15 
and 15.1 (Shi banov et alJ 12005*1 may be the 3rd and 
4th harmonics, respectively, which are fainter as the cy- 
clotron harmonic intensity decreases with its number n. 

Similar, albeit less significant, spectral features are 
likely seen in the photometric and spectral data of an- 
other middle-aged pulsar Geminga (Fig.l). The absence 
of strong nebular lines suggests that the features are also 
of the NS magnetospheric orig in. To explain the Geminga 
spectrum iMartin et alJ l)l998l) applied a toy model of an 
ion cyclotron absorption at B~10 n G in the inner mag- 
netosphere of the NS combined with the BB and PL 
components. At this field a low n ion cyclotron frequen- 
cies indeed fall in the optical range. However, it is likely 
that the nonthermal optical emission is generated in the 
upper magnetosphere where the magnetic field is by or- 
ders of magnitude weaker and any ion cyclotron absorp- 
tion of the respective optical continuum is negligible. In 
this case the electron/positron cyclotron absorption or 
scattering appears to be more plausible interpretation. 

The optical spectrum of the young Crab-pulsar is fea- 
tureless and has a different (positive) slope (Fig. 1). The 
ten times older Vela-pulsar has a flat and also featureless 
spectrum (Mignani, private communication). The spec- 



troscopy of PSR B0656+14 demonstrates that optical 
spectra of middle-aged pulsars can be distinct from those 
of younger ones by the presence of unusual spectral fea- 
tures or slope changes. To study this new spectral obser- 
vations of PSR B0656+14 in the NIR and 3000-5000A 
ranges are needed. A question, whether the observed dif- 
ference in pulsar spectra is simply caused by different 
pulsar geometry or by a change of physical conditions 
in the emission region with age, demands quantitative 
modeling physical processes in pulsar magnetospheres. 
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